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LIMITING SHAPE FOR DIRECTED PERCOLATION MODELS 

By James B. Martin 

CNRS and Universite Paris 7 

We consider directed first-passage and last-passage percolation on 
the nonnegative lattice Z+, d>2, with i.i.d. weights at the vertices. 
Under certain moment conditions on the common distribution of the 
weights, the limits g(x) = lim„._ 00 n~ 1 T( [nxj ) exist and are constant 
a.s. for x £ R+, where T(z) is the passage time from the origin to the 
vertex z £ Z+. We show that this shape function g is continuous 
on R + , in particular at the boundaries. In two dimensions, we give 
more precise asymptotics for the behavior of g near the boundaries; 
these asymptotics depend on the common weight distribution only 
through its mean and variance. In addition we discuss growth models 
which are naturally associated to the percolation processes, giving a 
shape theorem and illustrating various possible types of behavior with 
output from simulations. 

1. Introduction. We consider directed first-passage and last-passage per- 
colation models on the nonnegative lattice Zi, focusing in particular on 
behavior close to the boundaries of the orthant. 

With each node z £ Z+, associate the weight X(z). We assume that the 
hts {I(z),z£Zf} are i.i.d. according to some common distribution F 
on R; by allowing the weights to take negative as well as positive values we 
can consider first-passage and last-passage models simultaneously. 

A directed path in Z^ is a path each step of which consists of increasing a 
single coordinate by 1. Let T(z), the last-passage time to z, be the maximum 
weight of all directed paths from the origin to the point z, where the weight 
of a path is the sum of the weights of the nodes along the path. (See Section 
2 for precise definitions.) Natural objects of study are asymptotic quantities 
such as the function g : R defined by 

ET(Lnxj) 
S(x) =sup . 
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From superadditivity properties, we have that this supremum is in fact a 
limit, and that ?i~ 1 T([nxJ) — ► g(x) a.s. as n — > oo, for all x E R^. We call 
g the shape function, since it determines the limiting shape for the growth 
model associated to the percolation process; g(x) is also sometimes known 
as the time constant in direction x. 

Analogous first-passage percolation models on the undirected lattice are 
by now very well known — see, for example, Kesten [14, 15] or Durrett [7] 
for fundamental results. Recently, the directed last-passage model has also 
received much attention; in particular, the case where d = 2 and the weight 
distribution F is geometric or exponential. First, this is essentially the only 
nontrivial case (whether directed or undirected, first- or last-passage) where 
the form of the shape function g above is known; for exponential weights, it 
was first given by Rost [25]. But much more precise results are now avail- 
able; in particular Johansson [13] extended methods developed by Baik, 
Deift and Johansson [2] for the closely related model of the longest increas- 
ing subsequence of a random permutation, and showed that, for a > 0, the 
distribution of n~ 1 / 3 {T((n, [anj)) — ng((l,a))} converges as n — > oo to a 
nondegenerate limit (the Tracy-Widom distribution, which also arises as 
the limiting distribution for the size of the largest eigenvalue of a random 
matrix from the Gaussian unitary ensemble). 

Two-dimensional directed last-passage percolation problems with general 
weight distributions have also been studied in detail in the context of tandem 
queueing systems; see, for example, [1, 10, 19]. 

Our first observation is a condition on the weight distribution F under 
which the shape function g above is finite everywhere. The condition required 
on the positive tail is that / °°(1 — F(s)) 1 / d ds < oo; this follows quickly from 
analogous results for the related model of greedy lattice animals (introduced 
in [5] and [9]; the precise results we use are from [20]). We note that, as 
in the greedy lattice animals model, there is still a small gap between this 
sufficient condition and the best currently known necessary condition, which 
is that < oo (see [20] for a discussion). 

Our first main result is then that the shape function g is continuous on 
all of R^, including at the boundaries (in fact, continuity on the interior is 
immediate from a simple concavity property). We note that the question of 
continuity for the directed first-passage model was raised by Newman and 
Piza [21]; for the two-dimensional last-passage model it was resolved (in a 
queueing theory context) by Glynn and Whitt [10] for distributions with an 
exponential tail, and then by Baccelli, Borovkov and Mairesse [1] and Martin 
[19] under weaker moment conditions. Particular tools which we use to prove 
continuity at the boundaries in any dimension are a truncation which relies 
on a bound given in [20] for the growth rate in the greedy lattice animals 
model, and a concentration inequality derived from a result of Talagrand 
[27]. 
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In two dimensions we then give more precise information about the be- 
havior of g close to the boundary. For a distribution F with mean /x and 
variance <r 2 , we prove the asymptotic formula, as a — ► 0, 

g((l,a)) = n + 2<j\/a + o(y/a). 

In addition to the tools used in the proof of continuity, we use here a com- 
parison with a variant form of directed percolation analyzed by Seppalainen 
[26], and an estimate for the speed of convergence in the central limit the- 
orem from [24], in order to prove a universality property over all F for 
the asymptotics at the boundary Comparing the exact formula known for 
the case where F is the exponential distribution then yields the result. 
These asymptotics are linked to the Brownian directed percolation model — 
obtained, loosely speaking, by reversing the order of the limits n — > oo (in 
the definition of g) and a — > — which has been widely studied recently in 
various contexts; see, for example, [3, 11, 12, 23]. See also [22] for a sur- 
vey of the connections between these various directed percolation processes, 
random matrix theory and noncolliding particle systems. 

For various results on the dependence of the time constant on the weight 
distribution in the context of undirected first-passage percolation, see, for 
example, [16, 18, 28]. 

Just as in the case of undirected first-passage percolation, there are growth 
processes naturally associated to the directed percolation models considered 
here. In Section 5 we describe these and prove a shape theorem analogous to 
those given in [4] and [14]. We also discuss, with illustrations from simula- 
tions (see Figures 1-7), various possible behaviors of the growth processes, 
and the differences which exist between the directed and undirected cases 
and between the first-passage and last-passage cases. 

2. Notation and main results. We work with the d-dimensional nonneg- 
ative lattice Z+. For x £ 7L d + (and similarly Z d , M.+ and M. d ) we write Xi for 
the ith component of x; we use the norm ||x|| = \ x i\: an d write x < x' if 
x i < x i for i = 1, . . . , d. We write for the origin and 1 for the point all of 
whose coordinates equal 1, and for the point all of whose coordinates are 
except the ith which is 1 (so that x = Ylt=i x i e i)- 

With each point v of Z+, associate the weight X(y). We assume that 
the weights {X(v),v £ Z+} are i.i.d. random variables, with common dis- 
tribution function F, where F(s) = P(X < s) (as here, we sometimes write 
simply X to denote a generic random variable with distribution F). 

A directed path in is a path each step of which consists of increasing 
a single coordinate by 1. For zi, Z2 £ Z^, with z\ < Z2, let II[zi,Z2) be the 
set of directed paths from z\ to Z2. We identify a path with the set of points 
it contains; by convention we exclude the final point Z2 (but include the 



4 J. B. MARTIN 

initial point zi, unless z\ = Z2). Note that all paths in II[zi,Z2) have size 
(or "length") ||z2 — z\\\. 

For zi < Z2, define T(zi,Z2), the last-passage time from zi to Z2, by 

T(zi,z 2 ) = max V) *(v). 

In the case zi = 0, we write simply II[z) = II[0,z) and 

(2.1) T(z)=T(0,z) = max VX(v). 

7r£il z) 

1 ' VG7T 

We immediately have the following superadditivity property: if zi < z 2 < 
z 3 £ Zf, then 

(2.2) T(zi, z 2 ) + T(z 2 , z 3 ) < T(zi,z 3 ). 

Suppose E\X\ < 00. Then also E|T(z)| < 00 for all z G Z^. For x G R^, 
we now define 

(2.3) g(x) =sup-ET(L?ixJ) 

(which may be infinite). We sometimes write <7_f( x ) to emphasize the de- 
pendence on the distribution F; we also write g{x\, . . . for </(x) when 
x = (xi, . . . , xj) in order to avoid proliferation of brackets. 

The following basic properties of the function g are immediate from this 
definition and from the superadditivity in (2.2), using (a superadditive ver- 
sion of) Kingman's subadditive ergodic theorem: 

Proposition 2.1. Suppose E|X| < 00. 
(i) For oWxeK^, 
1, 



-T( [nxj ) — > <?(x) as n—> 00, a.s. and [«/ |<7(x)| < 00] in C\. 



(ii) 5 (ax) = a<?(x) /or a// a > 0, x G Rl. 

(iii) 5 is invariant under permutations of the coordinates. 

(iv) 5 (x) + g(y) < a(x + y) for all x, y G R d + . 

The following result gives conditions under which the function g is finite. 
Condition (2.5) is stronger than the condition that < 00 (which is 

known to be necessary for the finiteness) but weaker, for example, than the 
condition EA^(log + X) d ~ 1+e < 00. See [20] for details. 

Proposition 2.2. If 
(2.4) E|X|<oo 
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and 

(2.5) r \l - F{s)) 1/d ds < oo, 

Jo 

then \g{x)\ < oo for all x G R+. 
Our first main result is then: 

Theorem 2.3. Under conditions (2.4) and (2.5), the shape function g 
is continuous on all o/R+ (including at the boundaries). 

Proposition 2.2 and Theorem 2.3 are proved in Section 3. In Section 4 we 
prove the following theorem, which gives more precise asymptotics for g at 
the boundary in the case d = 2. 

Theorem 2.4. Let d = 2. Let the distribution F have mean fi and vari- 
ance a 2 , and satisfy 

(2.6) / (l-F(s)) 1/2 ds<oo 

Jo 

and 

r o 

(2.7) / F(s) 1 ' 2 ds < oo. 

J — oo 

Then as a j 0, 

(2.8) g(l,a)=n + 2ay^ + o(Va). 

Note that the framework effectively includes first-passage as well as last- 
passage percolation models, since the weights may take negative as well as 
positive values; replacing max by min and replacing the weights X(z) by 
—X(z) would simply change the sign of T and so of g. When considering 
associated growth models in Section 5, however, it is easier to consider first- 
passage and last-passage models separately. For completeness, we also state 
here the first-passage versions of the results above. Define the quantities 
{S(z),z G Z^_} and {/i(x),x G analogous to the last-passage quantities 
{T(z)} and {s(x)} defined at (2.1) and (2.3), by 

(2.9) S(z)= min VX(v), 

-en[z)^ 

(2.10) / l (x) = inf iE5(Lnxj). 

neN n 
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Corollary 2.5. (i) If E \X\ < oo and /^(^(s)) 1 ^ < oo, then \h(x)\ < 
oo for all x £ R^, and h is continuous on all q/R+ (including at the bound- 
aries). 

(ii) Let d = 2. If F has mean /i and variance a 2 and satisfies (2.6) and 
(2.7), then as a 10, 

h(l, a) = n — 2o\fa + o{\fa). 

Definitions, results and discussions for the growth models are given in 
Section 5. 



3. Continuity at the boundary. 



3.1. Case of bounded weights. We first prove the continuity result of 
Theorem 2.3 for the case where the weights are bounded. We will need the 
following concentration inequality, which follows easily from Theorem 8.1.1 
of [27]; see, for example, Lemma 5.1 of [20] for the argument. 

Lemma 3.1. Let Y%,i £ /, be a finite collection of independent random 
variables, such that 

F(\Yi\<L) = l 

for all i £ I . Let C be a set of subsets of I such that 

max I CI <R, 
cec ~ 



and let 



Then for any u > 0, 



Z = max > Yi . 

c&c 

i&C 



u 2 



n\Z - E Z\ > u) < exp (- + 64 

We apply the concentration inequality in the following lemma, which is 
the central part of the proof of the continuity of g. 



Lemma 3.2. Suppose F(\X\ < L) = 1 for some finite L. Let R > and 
e > 0. There exists 5 > such that i/x £ with ||x|| < R and Xj = (where 
1 < j < d), then 

\g(x + hej) - g(x)\ <e 

for all0<h<5. 
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Proof. Without loss of generality, let j = 1. Rephrased, the statement 
is that for x = (x2, 23, . . . , Xd) £ 

</((/i,x))->s((0,x)) as hiO, 

uniformly in {x: |x| < R}. 

So let h > and n G N. Any path from to the point ( \nh\ , [^xj ) contains 
exactly [nh\ steps which increase the first coordinate, so can be decomposed 
into a disjoint union of paths from (r, m r ) to (r, m r+ i), r = 0, 1, 2, . . . , [nh\ , 
where m r £ Zl" 1 for each r, and 



= m < mi < • • • < m Ln/l j +1 = [nxj . 
T(([nh\,[nx\)) 



(3.2) 



max 

m ,mi,...,m LnhJ+1 



(3.1) 

We have 

-\nh\ 

T((r,m r ),(r,m r+1 )) 

_r=0 

\nh\-l 
+ X(r,m r+1 ] 

r=0 

(The second term on the right-hand side appears because of the convention 
that a path from zi to Z2 does not include the "final point" Z2.) Here and 
below the m r range over values satisfying (3.1). The number of such choices 
for the m r is 

d 



n 

i=2 



[nxi\ + [nh\ 
[nh\ 



which, by Stirling's formula, is exp[n4>(h, x) + o(n)], where 
(j>(h,x)= Y ( h lQ g X% r - + x i lo s X 



2<i<d 
Xi>0 



h 



We now consider the expectation of the quantity within the maximum on 
the right-hand side of (3.2). For fixed {m r }, we have 

-\nh\ \nh\-l 

E ^ T((r,m r ),(r,m r+1 )) + X(r,m r+ i) 
_r=0 r=0 

= E ^T((0,m r ),(0, m r+l)) + L ra ^-JEX 
r=0 

(3.3) < ET((0,mo), (0,m^j +1 )) +nhL (by superadditivity) 
= ET((0, [nx]))+nhL 
<n[g((0,-x)) + hL] 
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< 



(by definition of g and superadditivity again). 

Still keeping {m r } fixed, note that the quantity inside the expectation on 
the left-hand side of (3.3) may be written as the maximum of the sum of 
various sets of weights X; each such set has size ||(Ln/iJ, L nx J)ll ^ n ll(^> x )ll 
and, by assumption, none of the weights has absolute value greater than L. 
So we can apply the concentration inequality in Lemma 3.1 to give 

-\nh\ \nh\-l 

T((r,m r ),(r,m r+1 ))+ ^ X(r, m r+ i) > n(g({0, x)) + hL + e) 

_r=0 r=0 

-\nh\ \nh\-l 

^r((r,m r ),(r,m r+1 ))+ X(r,m r+1 ) 

_r=0 r=0 

deviates from its expectation by more than ne 

( (ne) 2 
< exp ,,,, Nll n + 64 

f ne 2 
= eXP l-64||(/ l ,x)||L 2 +64 
Thus, taking the sum over all possible {m r }, 
F[T(([nh\ , LnxJ )) > n(g((0, x)) + hL + e)} 

( tie 2 

< exp[n*(M) + o(n)] ^{- 6m ^ L2 + 64 
This sums over n £ N to a finite amount whenever 



e > 8W||(/i,x)ll<KM). 

Since 

5 ((/i, x)) = hm -T(([nhJ , LnxJ )) a.s., 
Borel-Cantelli then gives 

g((h, x)) - g((0, x))<hL + 8W||(M)||<KM). 

The right-hand side tends to as h j 0, uniformly in ||x|| < R, as required. 

In the other direction, the superadditivity property in Proposition 2.1(iv) 
implies 

g((h,x))-g((0,x))>g((h,0)) 
= hKX 
> -hL, 

which again tends to as h j 0, uniformly over all x. □ 
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Lemma 3.3. Suppose W(\X\ < L) = 1 for some finite L. Then g is con- 
tinuous on all o/Mi. 

Proof. Let e > and y G W^. Suppose that y has exactly k nonzero 
coordinates. Without loss of generality, assume that in fact y« > for 1 < 
i < k and yi = for k + 1 < i < d. 

Define the function g k on by g k (u) = g(ui, . . . , u k , 0, . . . , 0) (appending 
d — k zeros to the end of u). 

Since g is concave on M.+ [by Proposition 2.1(iv)], g k is concave on Wl_, 
and so is continuous on the interior of R^L. 

Hence we can choose 5' > small enough that if \hi\ < 5' for 1 < i < k, 
then 



Choose any R > ||y|| + k5' . We now fix 5 > small enough that the con- 
clusion of Lemma 3.2 applies (for our chosen L, R and e), and also small 
enough that ||y || + kS' + (d — k)5 < R. 

Take any h G R d with ||h|| < min(<5',<5) and with (y + h) G M+. Then 
certainly \hi\ < S' for 1 < i < k, and also < hi < 5 for k + 1 < i < d. 

We are about to apply Lemma 3.2 (d — k) times, once for each of the 
coordinates of y which is 0. Specifically, for k + 1 < j < d, set x.^' = y + 
^i e «- Then for A; + 1 < j < d, we have ||x^) || < R (by choice of 5) and 

Xj = (since j/j = 0), so all the required conditions of Lemma 3.2 apply. 
Using also (3.4), we obtain 



g k (xi + hi,...,x k + h k ) 



g k {xi,...,x k )\ <e, 



and so 



(3.4) 







d 



+ J2 |5(xW + ^)-5(x^)| 



j=k+i 



<e + (d-k)e 
< (d + l)e. 
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Since e was arbitrary, we have that g is continuous at y, as desired. □ 

3.2. Extension to unbounded weight distribution. Define a lattice animal 
of size n to be a connected subset of Z, d of size n which includes the origin. 
Extend the i.i.d. array {X(z)} to all of 7L d , so that {X(z),z G is an i.i.d. 
array with common distribution F. Let A{n) be the set of lattice animals of 
size n, and define 

N(n)= max VX(z), 

the maximum weight of a lattice animal of size n. 

The results in the following proposition are taken from Theorems 1.1 and 2.3 
of [20]. (In fact the model in [20] covers only the case where the weights X(z) 
are nonnegative; but replacing weights whose value is by negative weights 
can only reduce the left-hand side of (3.5) or (3.6) and leaves the right-hand 
side unchanged. Alternatively, see [6] for a detailed treatment of the lattice 
animals model where the weights can take negative values.) 

Proposition 3.4. There exists c = c(d) < oo such that, for all F satis- 
fying (2.5): 

(i) for all n> 1, 

(3.5) EN(n) < cn [°° (1 - F(s)) 1/d ds; 

Jo 

(ii) with probability 1, 

(3.6) limsup^^ <c r°(l-F(s)) 1/d ds. 

n— >oo n Jo 

Now we can easily deduce the following lemma for the directed percolation 
model. Part (iii) implies Proposition 2.2. 

Lemma 3.5. There exists c = c(d) < oo such that, for all F satisfying 
(2.5): 

(i) for allzeZf, 

(3.7) ET(z)<c||z|| /°°(l-F(s)) 1/d ds; 

Jo 

(ii) with probability 1, 

(3.8) limsup- max T(z) < c f°° (1 - F(s)) 1/d ds; 

n^oo nz:\\z\\<n JO 
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(iii) for all x£R{ 



(3.9) 



x\\EX <g(x) < c||x 



r(l-F(s)) 1/d ds 
Jo 



Proof. First note that if z G and ||z|| = n, then any path ir G n[z) 
is a lattice animal of size n; thus T(z) < N(n), and parts (i) and (ii) follow 
immediately from Proposition 3.4. 

Putting z = [raxj in (i), dividing by n and then letting n — > oo gives the 
upper bound in (iii). 

For the lower bound in (iii), let z E Z^, and let tt be any path in n[z); 
then | vt" | = ||z||, and we have 



Again let z = [nx\ and let n — > oo to obtain the lower bound in (iii). □ 



We now introduce truncated versions of the weights {A^z)}. For L > 
and z G Z^, let X^(z) = max{min{X(z), L}, -L} [so that |A^ L )(z)| = 
mm(\X(z)\,L)}. 

Then let {T^fzJ.ze 1%) and {^(x),xe be defined just as {T(z)} 
and {g(x)}, but with the quantities {X(z)} replaced by the truncated ver- 
sions {X^ L \z)}. 

Lemma 3.6. Suppose that (2.4) and (2.5) hold. Then for any xEffi+, 




= llzllEX 



!J 




— x 




(3.10) 




where c is as in Lemma 3.5. Thus, for any R > 



(3.11) sup | 9 (x)- 9 (t) (x)h0 



as L 



CO. 



x£lH : ||x||<R 



Proof. 



Note that 



[L - X(z)} + < X(z) - I«(z) < [X(z) - L} + . 
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We consider first the positive tail. Let x G . Then 
g(x) -# {L) (x) = lim -ET(|_nxJ) - lim -E ( [nx] ) 



(3.12) 



ri — *-oo 77, 

lim — E 

n— >oo 

1 



n — >oc ji 

sup ^X(v)- sup J]X (i) (v) 

_7ren[0,LnxJ) vg7r 7ren[0,LnxJ) ve7r 



< lim — E sup 

™^°°n we n[o,Lnxj) 



£a»-£a^)(v) 



.VG7T 



V€7T 



lim — E sup 

,H0On -Gn[0,Lnxj 



£[X(v)-L]- 



< C X 



(l-FOO) 1 /^; 



the last inequality follows from Lemma 3.5, since the variables {LY(v) — L] + , 
v G Z^} are i.i.d. with distribution function F^ >L \ where F( >L \s) = for 
s < L and F^ >L \s) = F(s) for s > L. This gives the lower bound in (3.10). 

For the negative tail, let z G Z+, and let 7r* G II[0,z) be the maximizing 
path for T^(z). (If there are several maximizing paths, choose, say, the one 
that is first in the lexicographic order.) 

Now for any v G Z+, P(v G 7r*|X(v) < s) is a nondecreasing function of s. 
[This follows from a simple coupling, since {X(V), v' G Zl} are independent 
and, for a fixed realization of the other weig hts {X(v')y €Z%y t^v}, the 
function /{v G 7r*} is a nondecreasing function of X(v).] 

Hence in particular P(v G ir*\X(v) < — L) < P(v G vr*), and so by simple 
manipulation of conditional probabilities, 

P(X(v) < -L|v G vr*) < P(X(v) < —L). 

Furthermore, the event {v G tt*} depends only on {max{X(v'), — L}, v' G 
Z^f_}; thus, conditional on {X(v) < —L}, X(v) is independent of {v G it*}. 
Hence 

E([-L-X(v)] + |vG7T*) 

= E ([— L - X(v)]+|X(v) < -L)P(X(v) < -L\v G vr*) 
< E ([—L - X(y)]+\X(v) < -L)F(X(v) < -L) 
= E([-L-X(v)]+) 



-L 



F(s)ds. 
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Now 



ET(z)=E max VX(v) 

^Wveir 

>E max V(xW(v)-[-L-X(v)] + ) 

L ' VG7T 

>E ^#(v)-E £[-L-*(v)]+ 

VG7T* VG7T* 

= ETW(z)- £ P(vGvr*)E([-L-X(v)] + |vG7r*) 
>ET (L) (z)-/ ^ P(vevr*) 

= ET^(z) - / F(s)dsE|7r*| 

»/ — oo 

r—L 

= ET^(z) - ||z|| / F(s)ds, 

J — oo 



since tt* £ H [z) and all paths in II [z) have length ||z||. 

Letting z = [^xj and taking n — > oo then gives the first inequality in 

(3.10) , as required. The convergence in (3.11) then follows since, under (2.4) 
and (2.5), both the integrals in (3.10) tend to as L —> oo. □ 

It is now immediate to extend the continuity property proved in Lemma 
3.3 to the case of unbounded weights, and so complete the proof of Theorem 
2.3. 

Lemma 3.3 shows that the functions g( L ) are continuous for each L, and 

(3.11) shows that, under the conditions (2.4) and (2.5), — ► g as L — > oo, 
uniformly on any compact subset of Hence g itself is continuous, as 
required. 

4. Asymptotics at the boundary for d = 2. In this section we prove The- 
orem 2.4. 

We first obtain an estimate on the growth rate in the case where F is a 
Bernoulli distribution. This is done using a comparison with an alternative 
percolation model in which the Bernoulli distribution is an exactly solvable 
case. We write Ber(j>) for the Bernoulli distribution with parameter p, with 
p(X = l) = l-P(X = 0)=p. 

Lemma 4.1. For all a > 0, pG [0,1], 



ffBer(p) (1, a) < (1 + a)p + 2y/a\/l + aVp(l-p). 
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Proof. Recall that 




z 0> z l i • • • i z m+rt— 1 



such that zq = 0, such that, for all 1 < i < m + n — 1, Zj — z^_i = e,- for 
some j G {1, . . . , c£}, and such that also (m,n) — z m+n „i = for some j G 



Define an alternative set of increasing paths U[m,n) to be those paths of 
the form 



Define the function ijj : i— ► by tp(x, y) = (x + y, y). 

Since VK 2 ) 7^ VK 2 ') whenever z 7^ z', and since {X(z),z G Z?_} are i.i.d., 



(4.1) KT(m,n) = E max VX(z)=E max VJT(^(z)). 

7ren[m,n) ze7r 7r61I[m,n) 

For a path 7r = zo, • ■ • , z r _i, write V^ 71 ") f° r the path ^(zo), . . . , ^(z r _i). 
From the definitions jaf the path sets above, one can obtain that if tt G 
U[m,n), then i/j(tt) G n[m + n,n). (We do not require that all paths in IT[m + 
re, n) can be written as VK 71 ") m this way.) Put another way: let tp(H[m, n)) be 
the set {fr : tt = for some 7r G n[m, n)}; then ^(n[m, n)) C n[m + n, re). 

Continuing from (4.1), 



{!,... 



d}. 



(0,yo),(l,yi),...,(m- l,y m -i) 
where < yo < y\ < ■ ■ ■ < y m ~i < re, and define 

T(m,n) = max V^.X'(z). 

7ren[m,n) zG ^ 



we have 



ET(m, 




< E jnax ^ X(y) 



n£H[m+n,n) v g^- 



= ET(m + n, re). 



Then 



y(l,a) 



lim — ET(n, I an 



(4.2) 



< 



liminf — ET( [1 + ajre, |_cmj). 
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Seppalainen [26] analyzes directed percolation based on the path sets II, 
and in particular obtains that, for the case of Bernoulli weights, 

lim -E Bcr / p) f ([aiu], [a 2 n\) 

n— >oo n 

a i 



p(ai - a 2 ) + 2^a 1 a 2 \/p(l - p), itp< 

Oi\ + «2 

Oil, it P> • 

A calculation then shows that for all p, 

\un o ^E Bcr ^f([a 1 n\, [a 2 n\) < a.\p + 2y/aia 2 Vp(l -p). 

Substituting into (4.2) with a\ = 1 + a, a 2 = a gives the required result. □ 

Lemma 4.2. Let Fx and F 2 be distributions with means p\ and fi 2 and 
satisfying (2.6) and (2.7). Then for all a > , 

\g Fl {l,a) -g F2 (l,a) - (1 + a) (fix - // 2 )| 

< 2*/a(l + a) r \F 1 (s)-F 2 (s)\ 1/2 ds. 

J — oo 

Proof. Let {U(z),z G be i.i.d. uniform on [0,1], and for i = 1,2, 
let Xi(z) = Fr 1 (U(z)), where F _1 (-u) = sup{x : F(x) < u}. 

Then for i = 1,2, {JQ(z),z G are i.i.d. with distribution F^, and for 

any x, 

P(Xi(z) > x,X 2 (z) < x) = [F 2 (x) - 

and 

P(X 2 (z) > x,Xi(z) < x) = [Fx(x) - F 2 (x)}+. 

Now 

g Fl (l,a) -g F2 (l,a) 

= lim — E max y^Xi(z)— lim E max X 2 (z) 

n»Jl 7ren(n,|cm|) V n->oo TrelKn, I an I ) ^ 

J / z^z-tt J ' ZG7T 



< lim -E max Vp^fz) - X 2 (z)) 

= lim -E max / V[J(Xi(z) > x,X 2 (z) < x) 

rwoon Trenfn.LanJ) J -oo ^ 

- J(Xi(z) <x,X 2 (z) >x)]dx 
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l r°° 

< lim — E / max VfJ(Xi(z) > x,X 2 (z) < x) 

-n^oon 7_oo7ren(n,Lanj)^ L W ~ 

-/(Xi(z) <x,X 2 (z) >x)]dx 

/oo J 
lim — E max V fI(Xi(z) > x, Jf 2 (z) < x) 
-oo™'! Ten(n,Lanj)^ L 

- I(Xi(z) < x,X 2 (z) > x)) dx 
(by Fubini's theorem and bounded convergence) 



/oo if 
lim -^E max V llfXUz) > x,X 2 (z) < x)] 
-oo^n) Ten(n,Lanj)^ L 

(^[l-/(X 1 (z)<x,X 2 (z)>x)] 
\ze7r 



+ E max 

7rGiI(n, [cmjt) 



J2i)\dx 

ze-rr / ) 



{9Bcv([F 2 {x)-F 1 (x)] + )( 1 i a ) + 9Ber(l-[F 1 (x)-F 2 (x)] + )(' i -i a ) ~ (1 + Ot) dx} 

oo 

/oo 
{(1 + a)([F 2 (x) - + 1 - [F(x) - F 2 (x)]+) 

-oo 

+ 2^VTT^([F 2 (x) - F 1 (x)] 1 / 2 + [Fi(x) - F 2 (x)]V 2 ) 

- (l + a)}dx 

(by Lemma 4.1) 

/oo 
-oo 

/oo 
(^(x)-^^)! 1 / 2 
-oo 

/oo 
|Fi(x) -F 2 (x)\ l ' 2 dx. 
-oo 

Similarly, 

5F 2 (l,a) - g Fl (l,a) 

/oo 
|Fx(x) -F 2 (x)| 1/2 dx. 
-oo 

Together these give the desired result. □ 
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Lemma 4.3. Let F satisfy (2.6) and (2.7), and let e > 0. Then there is a 
distribution F with bounded support which has the same mean and variance 
as F , and which satisfies 



\F{s)-F{s)\ 1 ' 2 ds 



< e. 



Proof. Let X have distribution F. For brevity we cover only the case 
where F(X > 0) = 1; the negative tail can be truncated in an analogous way. 

Take any t > 0. If ¥(X > t) = 0, then F itself has bounded support and 
we take F = F. Otherwise, let m = E (X\X > t) and w = E (X 2 \X > t), and 
choose p, u to satisfy 

(4.3) (1 — p)t + pu = m, 

(4.4) (1 -p)t 2 +pu 2 = w. 
The solution is 

(m-t) 2 

P 



(m — t) 2 + w — m 2 ' 
m — t 



u = t + 

P 

note that < p < 1 and u> t, since m>t and w > m 2 . 
Then let 

( F(x), if0<x<t, 
(4.5) F(x) = <l-p[l-F(t)}, ift<x<u, 

[ 1, if x > u. 

Now F has bounded support since F(u) = 1. For the mean and variance, 
let X have distribution F, and use (4.3) and (4.5) to give 

EI = E (X;X < t) + (1 - p)[l - F(t)]t + p[l - F(t)]u 

= E (X; X < t) + P(X > t)[(l - p)t + pu] 

= EX; 

similarly use (4.4) to give EX 2 = EX 2 . 
For the final part we have 

/ \F(s)-F(s)\ 1/2 ds< / [l-F(s)] 1/2 ds+ / [1 - F(s)} 1/2 ds 
Jo Jt Jt 



= J [1-F{sp 2 ds + J [p{l-F{t))] 1/2 ds 

< J [l-F(s)] 1 / 2 ds+ \pu 2 P(X > t)] 1 / 2 

POO 

< J [l-F(s)] 1/2 ds + E(X 2 ;X>t) 1 / 2 . 
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By assumption, f£°[l — F(s)] 1 ^ 2 ds < oo, and this implies that EX 2 < oo 
also; hence by choosing t large enough we can make the right-hand side as 
small as desired. □ 



Lemma 4.4. Let F be a distribution with bounded support, and, for k G 
N, let F^ be the distribution of X\ + X2 + • • • + X^, where {Xi} are i.i.d. 
~ F. Let r :R + 1— > N be any function satisfying r(a) — > 00 and r(a)y / a — ► 
as a J, 0. TTien 

of(1, a) - Z7Z\9F(r) (!) ar(a)) = 0. 



lim — p= 

«|0 



r(a)' 



Proof. For (x,y) G Z+ and r G N, let B^ r \x,y) be the set {(rx + 
«,y),i = 0,l,...,r- 1}. 

The sets B^ r '(z) partition Z+; essentially we have grouped the sites of Z^j_ 
into "blocks" of length r and height 1. We will compare our original model 
with one where each of these blocks functions as a single site, whose weight 
is the sum of the original sites contained in the block. 

Given tt G H(nr,m), we can find tt G TL(n,m) such that 



(4.6) 



[J B (r) (z)A 



< mr, 



where A denotes the symmetric difference. [For instance, let tt = {z : tt n 
B^~\z) 7^ 0}.] Similarly, given tt G Tl(n,m), we can find tt G H(nr,m) such 
that (4.6) is satisfied. 

Suppose that {X(z),z G Z 2 ^} are i.i.d. with distribution F. Define 

x^{z)= Yl X ^)- 

v£flW(z) 

Then {X» (z), z G Z 2 } are i.i.d. with distribution F^ r ' . Let if be such that 
P(|X| > K) = 0. Then by the properties (in both directions) noted at (4.6) 
and after, 



max y^X(z)— max Y"' X(z) 

n£U(nr,m) _ 7r£in(n,m) 



< mrK, 



so that 



1 11 

— E pT(nr, m) E F ( r )T(n, m) 

nr r n 



< — mrK. 

nr 



Putting m = [anr\ and letting n — > 00 gives 

gF(l,a) - -g F (r){l,ar) 
r 



< arK. 
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If r is a function of a such that rJa — > as a j 0, then the right-hand side 
is o(\/a) as a | 0, as desired. □ 

Lemma 4.5. Let F be a distribution with bounded support, with mean 
fip and variance op. Let F^ and r be as in Lemma 4.4. Then 

g$(l,ar(a)) 



lim —= 

o|0 Jet 



1 



rip) 



g FM (1, ar(a)) - fi F - <r F - 



o. 



where is the standard normal distribution. 



Proof. Theorem 5.16 of [24] gives a bound on the rate of convergence 
in the central limit theorem, for distributions F which have a finite third 
moment. Here F has bounded support and so certainly finite third moment; 
we obtain that there exists C = C{F) such that 

\F^{x) - $(x)| < Cr^ 2 (l + \x\)~ 3 

for all r G N, x £ R, where F^ r > is the distribution of (X\ + • • • + X r — 
r/j, F )/(o~ F Jr). (Note that F*) is simply the distribution F^ normalized to 
have mean and variance 1.) 

Now combine this estimate with Lemma 4.2; for any r £ N, 

11 . 0$(l,ar) 

-g F (r)(l,ar) - fi F - v F -= — 

r x/r 



1 

7 a 
op 



o- F 5$(l,ar) 
-j=g~ {r) {l,ar) - o-p -=— 



gp {r) (l,ar) - g$(l,ar)\ 



< -^L2x/a7VT+a~r~ f°° iCr^n + |^|) — 3 1 1/2 dx 
Jar J-oo 

a/1 + ctr 



C'a F - 



,1/4 



where C is some constant independent of r and a. If r is a function of a 
such that r\fa — ► and r — > oo as a [ 0, then the right-hand side tends to 
as a | 0, as required. □ 

The following lemma is the universality result which we need: 

Lemma 4.6. Let F be a distribution with mean fi F and variance a F , 
and satisfying (2.6) and (2.7), and let the function r satisfy the conditions 
of Lemma 4.4. Then 



(4.7) 



lim —= 

aj.0 \Ja 



g F (l,a) — /j, f — o~p 



g<s>(l,ar(a)) 



x/rja) 



0. 
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Proof. Let e > 0. Using Lemma 4.3, choose a distribution F with 
bounded support, with the same mean and variance as F, and with 

1/2 

-oo 



(4.8) 
Then 



limsup — -j= 

aj.0 V" 



\F(x)-F{x)\ L/z dx<e/2. 

g$(l,ar(a)) 



g F (l, a) - [i F - <Jf- 



1 

< hm sup —= 



«) - Hf-<?f 



x/rta) 

5$(l,ar(a)) 



1 



+ limsup -^=1^(1, a) - c^(l, a) 



a|0 



The first term is by combining Lemmas 4.4 and 4.5 and using the fact that 
F has the same mean and variance as F; the second term is < e by Lemma 
4.2 and (4.8). This works for any e > 0, so the desired result follows. □ 

Finally, we compare with an exactly solvable case to yield the asymptotic 
behavior for all F: 

Proof of Theorem 2.4. Choose any r that satisfies the conditions of 
Lemma 4.4, for example, r(a) = L^" 1 ^ 4 ] ■ 

When F is the exponential distribution with mean 1 (and so also vari- 
ance 1), we have the exact formula a) = 1 + 2-^a + a (see, e.g., [25]). 
Substituting into (4.7) gives 

5$(l,ar(a)) 



(4.9) 



lim —= 

aj.0 y/a 



0. 



Now take any F satisfying (2.6) and (2.7). Combining (4.7) and (4.9) 



lim— =|<7i?(l, a) — llf — 2aF\fa\ = 0, 
«J.o J a 



gives 



as required for (2.8). □ 
5. Growth models. 

5.1. Definitions and statement of shape theorem. Recall that first-passage 
quantities S and h, analogous to T and g, were defined at (2.9) and (2.10). 
Define B(t), the last-passage shape at time t, by 



B(t) = {xeM d + :T([x\) <t}, 
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and define C(t), the first-passage shape at time t, by 

C(t)={xeR d + :S([x\)<t}. 

Both B(t) and C(t) are increasing in t, in the sense that for < t\ < <2, 
B(h)QB(t 2 ) and C(ii) C C(t 2 ). 

We further define subsets B and C of M.i by 

B = {x:(/(x)<l}, 

C = {x:/i(x) < 1}. 

B is concave [by Proposition 2.1(iv)]; similarly, C is convex. B and C are 
asymptotic shapes for the processes {B(t)} and {C(t)} in the sense of the 
following theorem, which is analogous to well-known results for undirected 
first-passage percolation models (see, e.g., [4, 14]). We give the proof in 
Section 5.3. 

Theorem 5.1. Suppose that the weight distribution F satisfies 

(5.1) f° F(s) l ^ d ds < oo 
and 

POO 

(5.2) / (l-F(s)) 1/d ds<oo. 

Jo 

(i) Last-passage shape theorem. 
If E pX > 0, then for any e > 0, 

(5.3) (l-e)BC^C(l + £ )B 

for all sufficiently large t, with probability 1. 

(ii) First-passage shape theorem. If 

, h(x) 
inf -f-f > 0, 

xGMl\{0} ||x|| 



then for any e > 0, 



;i-£)CC^C(l + e)C 



for all sufficiently large t, with probability 1. 

Remark 5.1. (i) Note that by a subadditivity property for h, analo- 
gous to the superadditivity property for g in Proposition 2.1(iv), we have 
that for all x G R 1 ^, h(x) > ^-h(l, !,...,!); thus the condition in part (ii) 
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is equivalent to the condition that h(l, 1, . . . , 1) > 0. If the weights are non- 
negative, then this is implied by the condition that F(0) <Pc, where 
is the critical value for directed percolation in d dimensions; this follows, for 
example, from the same arguments as the property, noted by Kesten and 
Su [17], that (in their case for undirected percolation) the critical points for 
percolation and for "1-percolation" coincide. 

(ii) In fact, the moment conditions in Theorem 5.1 are stronger than 
necessary; for the last-passage case one can replace (5.1) by the condition 
that E l-X—l < oo, and for the first-passage case one can replace (5.2) by the 
condition that EX < oo. 

In particular, combining with the previous remark, for the first-passage 
model with nonnegative weights it suffices for the limiting shape result that 
EX<oo and F(0) <p { c d) . 

To prove the theorem under these weaker conditions, one can follow an 
approach similar to that used by Cox and Durrett [4] for the undirected first- 
passage case, making use of the fact that E X < oo E min(Xi, . . . , Xd) d < 
oo, where X%, . . . , are i.i.d. copies of X. However, the proof needs a rather 
lengthy enumeration of cases and a description of a variety of different sets 
of "alternative paths" ; in addition, the fundamental ideas are already in [4] , 
so we do not include it here. 

Under the stronger conditions in Theorem 5.1 (already almost optimal for 
the last-passage case with nonnegative weights) , bounds of the sort given in 
Lemma 3.5 are available, and the proof is simpler. 

5.2. Discussion of the growth models. In this section we describe and 
illustrate various possible types of behavior for the first-passage and last- 
passage growth processes B(t) and C(t). 

We will concentrate on the case where the weights {X(z)} are nonnega- 
tive. Then for all t > 0, B{t) and C(t) are connected; in addition, B(t) is a 
decreasing subset of , although C (t) does not generally have a similarly 
simple property. 

If the weight distribution F is exponential (resp. geometric), then the pro- 
cesses {B(t),t > 0} and {C(t),t > 0} are Markov in continuous (resp. dis- 
crete) time; in the two-dimensional last-passage case this yields the growth 
model considered in, for example, Rost [25] and Johansson [13]. Simulations 
of B(t) and C(t) in two dimensions with exponentially distributed weights 
are given in Figures 1 and 2, and the three-dimensional last-passage case is 
shown in Figure 7. 

In fact B(t) and C(t) are also Markov in discrete time when the weights 
are Bernoulli (taking values and 1). 

We now discuss various ways in which the shape result in Theorem 5.1 
may fail. 
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First, the last-passage case. Note that if g(x) = oo for some x in the inte- 
rior of Ri, then (by a simple superadditivity argument) g = oo throughout 
the interior. A sufficient condition for this to occur is that EX d = oo. Then 
the growth of B(t) in any interior direction is sublinear in t; on a linear scale, 
the asymptotic shape collapses into the boundary (or even to the origin alone 
if EX = oo). An example is illustrated in Figure 3 for a distribution with 
finite mean but infinite variance, with d = 2. 

The undirected first-passage case on Z d is comparable. Here again there 
are just two possibilities: either the shape function is everywhere (in which 
case the asymptotic shape is essentially the whole of R d ), or the shape 
function is nonzero everywhere. 

For the directed first-passage case, in contrast, the different behaviors 
can co-exist. If < ¥(X = 0) < 1, then h = for some cone around the 
direction 1 (this is the cone in which "oriented percolation of sites with 
weight 0" occurs), but h > elsewhere. In this case there will a.s. be some 
infinite path starting at the origin which has finite total weight, and the set 
C(t) will have infinite size at some finite time. The shape C = {x: h(x) < 1} 
is noncompact, but not equal to the whole of Ml, and the shape theorem 
does not apply as given. See Figure 4 for a simulation of such 

Note that if again X attains its minimum value with probability more than 
pi d \ but this minimum is now greater than 0, then /i(x)/||x|| is constant on 
some cone around the direction 1, and the boundary of C has a flat section — 



600 1 




Fig. 1. Simulation of the last-passage 
process for d = 2 and F exponential with 
mean 1. The sets B(t) are shown for 
t=150 (darkest), 300,450,600 (lightest). 
Here the asymptotic shape B is known to 

fe{Mei^+Vj<i}. 




Fig. 2. First-passage process C(t), 
t = 50,100,150,200 for d = 2 and 
F exponential with mean 1. 
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see Figure 5. This is the same phenomenon observed by Durrett and Liggett 
[8] for the undirected case. 

If EX = oo, then h(ei) = oo, and the shape theorem cannot hold as stated, 
but an amended version for a cone excluding the boundaries may hold. 
However, if (only slightly more strongly) E min(Xi, . . . , Xd) d = oo, then the 
result fails more fundamentally. Still the limit /i(x) = lim rwoo n _1 S , (L?i-xJ) 
exists and is finite and constant a.s. for any x in the interior of Ri, and 
one can define the asymptotic shape C as before. However, it is no longer 
the case that the convergence is a.s. uniform on compact subsets of in 
effect, the "holes" in the shape C(t) persist, as seen in Figure 6. The same 
sort of behavior would occur for the undirected first-passage model when 
E min(Xi, . . . ,X2d) d = oo; see, for example, related discussions in [4]. 

5.3. Proof of Theorem 5.1. Note first that for part (i) of Theorem 5.1, 
the condition E X > implies that inf xgR d j^- > [since, by superad- 
ditivity, g(x) > \\x\\g(l, 0, . . . , 0) = ||x||EX]. 

Note also that, by replacing the weights {X(z)} by {— X(z)}, part (ii) 
can be rewritten as follows in terms of last-passage rather than first-passage 
quantities: 




Fig. 3. B(t), t= 150,300,450,600, with Fig. 4. C(t), t = 18,36,54, with d = 2 

d — 2 and the Pareto distribution and X = w.p. 0.645 and X = (0.355) -1 

F(x) — min(0, 1 — (3x) _3,/2 ), which has w.p. 0.355 (showing only the intersection 

mean 1 but infinite variance. The asyrnp- with the box [0,350] 2 ,). Here /i(x) = for 

totic shape B consists only of two lines, some (though not all) x; thus the asymp- 

between the origin and (0, 1) and between totic shape is noncompact (but not the 

the origin and (1, 0). whole o/R+j, and a.s. the set C(t) will be 

infinite for some t. 
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50 100 150 200 250 300 




20 40 60 80 100 120 140 



Fig. 5. C(t), t = 50, 100, 150, 200, with Fig. 6. C(400) ford = 2 and the Pareto 
d—2 and X = 0.5 w.p. 0.8, X — 3 w.p. distribution F(x) — min(0, 1 — x~ 3 ^ 4 ), 
0.2. which has infinite mean. 
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Fig. 7. B(30) for d — 3 and F exponential with mean 1. 



If 



sup ^<0, 

xgM^\{0} ll x ll 



26 



J. B. MARTIN 



then for any e > 0, 

(5.4) (l-£)rc^C(l+ £ )r 

for all sufficiently large t, with probability 1, where 

B-{t) = {xeR d + :T{lx\)>-t} 

and 

£T = {x: 5 (x)>-l}. 

By Theorem 2.3, we know that |<?(z)| < oo for all z; then (5.3) and (5.4) 
are immediately implied by the following property: for any e > 0, there are 
a.s. only finitely many z G Z^ such that \T(z) — g(z)\ > e|</(z)|. 

Since we assume in both cases that |g(z)|/||z|| is bounded away from 0, 
we have |<?(z)| > ||z|| inf z / ^(z'^/Hz'H, so in fact it is enough to show that for 
any e > 0, there are a.s. only finitely many z such that |T(z) — g(z)\ > s\\z\\. 

This follows immediately from Lemmas 5.2-5.5, which we state immedi- 
ately and then prove in turn. 

Lemma 5.2. Suppose F satisfies F{s) 1 t d d8 < oo and / °°(1 - F(s)) 1 ^ ds < 
oo, and let e > 0. If L is large enough, then with probability 1, 

|T(z) -T( L )(z)| <e||z|| 

for all except finitely many z G Zi . 

Lemma 5.3. Let e > and L > 0. With probability 1, 
|TW(z)-ET( x )(z)| < e ||z||, 
/or aZZ except finitely many z G Zl . 

Lei e > and L > 0. T/ien /or aZZ except finitely many 

\ET^ L \z)-g^ L \z)\ <e||z||. 

Lemma 5.5. Suppose F satisfies f® 00 F{s) l l d ds < oo and / °°(l-F(s)) 1 / d ds < 
oo, and let e > 0. 7/ L is sufficiently large, then for all z G Zi, 

|^ L )(z)- 5 (z)| <e||z||. 



Lemma 5.4. 
zgZ^., 
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PROOF of Lemma 5.2. Let L be large enough that c — F(s)) l l d ds < 
e/2 and c JZ OQ F(s) 1 ^ d ds < e/2, where c is the constant in Lemma 3.5. 
Let z £ Z^; for some 7r* G n[z), we have 



T(z)-rW(z)= £[X(v)-A^)(v)], 

Vg7T* 

so that 

|T(z) - #)(z)| < 2 [*(v) ~L]++ E [" L " ^( v )]+ 

VG7T* VG7T* 

(5.5) 

< max V F (L) (v)+ max V W (I,) (v), 
~7ren[ z )^ 7ren[ z )^ 

1 ' VG7T L ' Vg7T 

where we define V^ L \v) = [X(y) - L]+ and W^ L \v) = [-L - X(y)]+. 

Note that {W L )(v),v G Z^_} are i.i.d. with common distribution Fy\ 
where Fy 1 ^ (x) = for x < and (x) = F(L + x) for x > 0. 

Similarly, {#)(v),v£Z|} are i.i.d. with common distribution F^\ 

where F^\x) = for x < and F^\x) = 1 — F(—L — x) for x > 0. 

From Lemma 3.5 (applied to {V( L \v)} and Fy 1 ^ rather than to {X(v)} 
and F), we then have that, with probability 1, 



limsup— max max S~] V^iy) < c [ (1 — Fy L \s)) l ^ d ds 

n— >oo n z: ||z||<n7ren[z) Jo 

:J™(l-F(s)f d 



= c 
<e/2. 

In particular, there are a.s. only finitely many z G Z+ such that 

max V l/ (L) (v) > -llzll. 

1 ' VG7T 

Applying Lemma 3.5 to {W^(v)} and F^ in the same way, one obtains 
that there are a.s. only finitely many z G Zi such that 

max V W {L) (v) > -llzll. 

1 ' VG7T 

Thus from (5.5), there are a.s. only finitely many z G Z+ such that 

|r(z)-r( L )(z)|> £ ||z||, 

as required. □ 
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Proof of Lemma 5.3. All the paths in n[z) have length ||z||, and the 
weights X( L '(y) have absolute value no greater than L. Hence we may apply 
the concentration inequality in Lemma 3.1 to give 

P(|TW( Z ) _ ET( L )(z)| > e ||z||) < exp(-|M^ + 64 

=exp (-S^ +64 )- 

For re G Z+, there are certainly no more than (n + l) d points z such that 
||z|| = n; thus 

£ P(|r( i )(z)-ET( i )(z)|>e||z||)< £ (^ + l)"exp(-^ + 64 

< oo, 

and Borel-Cantelli yields the desired result. □ 

Proof of Lemma 5.4. From the definition of g and so of g^ L \ we have 
ETW(z) < g {L) (z) for all z, so we need to show that ET^fz) > # (L) (z) - 
e||z||, except for finitely many z. 

Fix a > 0. The distribution has bounded support, so certainly sat- 
isfies (2.4) and (2.5); thus by Theorem 2.3, g^ is continuous on Ri, and 
hence is uniformly continuous on the compact subset {x G : ||x|| < 2d}. 

So choose u < min(l,a) such that whenever ||x|| < d and ||x — x'|| < ud, 
then |^ (L) (x) -5 (L) (x')| < a. 

Now let C be the set 



ur, r£ < 0, 1, . . . , 



1 

u, 

C is a finite subset of M+, and for every y G C, we have [by Proposition 
2.1(i)], 

] (y) asn^oo. 

re 

Hence there is iV = -/V(a) such that, for all n > N and all y G C, 

ET^([ny\)>n(g^(y)-a). 
Let z satisfy max Zi > iV. Define 

1 z 



u max z,- 
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then y G C, with (maxzj)y < z, with ||y|| < d and with 



<ud< ad. 



max z,- 



Using first superadditivity, then the fact that all the weights {X^ L \z)} 
are no smaller than —L, then the continuity bounds above, we obtain 

ETW(z) >ET^(L(maxzi)yJ) +ET ( - L \z- [( maxz; y ) 

> ET( L )(L(maxZ J )yJ) - L\\z - [(maxZ^yJ || 

> (max^)(5( (L) (y) - a) - L(\\z - (max^)y|| + d) 



# (L) (z) - (maxzj) 



+ a + L 



\ max Z{ 
z 



9 iL Hy) 



max z,- 



+ 



Ld 



max z,- 



> (z) — (maxzj)|a + a + Lad + 



Ld 



max z,- 



Hence if a < 2e 1 (2 + Ld), then for all z with maxzj > max(iV(a), 2Ld/e), 
we have 

(maxzj)e 
>g {L) (z)-e\\z\\, 

as required. □ 

Proof of Lemma 5.5. Under the moment conditions on F, the result 
follows immediately from Lemma 3.6. 

This completes the proof of Theorem 5.1. □ 
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